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=3 ■ ABSTRACT 

S: 

^ ■ We confront predictions of inflationary scenarios with the WMAP data, in 

combination with complementary small-scale CMB measurements and large- 

^ ■ scale structure data. The WMAP detection of a large-angle anti-correlation 

^ . in the temperature-polarization cross-power spectrum is the signature of adi- 

CSl ■ abatic superhorizon fluctuations at the time of decoupling. The WMAP data 

^ . are described by pure adiabatic fluctuations: we place an upper limit on a cor- 

. related CDM isocurvature component. Using WMAP constraints on the shape 

^ . of the scalar power spectrum and the amplitude of gravity waves, we explore 

Q_|. the parameter space of inflationary models that is consistent with the data. 

Q . We place limits on inflationary models; for example, a minimally-coupled X(j)^ 

^ , is disfavored at more than 3-0" using WMAP data in combination with smaller 

. scale CMB and large scale structure survey data. The limits on the primordial 

> ! parameters using WMAP data alone are: Usih = 0.002 Mpc~^) = 1.20toJ^, 
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dus/dhik^ -0.0771^:^^^, A{ko = 0.002 Mpc^^) = O.llto.n (68% CL), and 
r{ko = 0.002 Mpc-^) < 1.28 (95% CL). 

Subject headings: cosmic microwave background — cosmology: observations — 
early universe 

1. INTRODUCTION 

An epoch of accelerated expansion in the early universe, inflation, dynamically resolves 
cosmological puzzles such as homogeneity, isotropy, and flatness of the universe (Guth 1981; 
Linde 1982; Albrecht & Steinhardt 1982; Sato 1981), and generates superhorizon fluctuations 
without appealing to fine-tuned initial setups (Mukhanov & Chibisov 1981; Hawking 1982; 
Guth & Pi 1982; Starobinsky 1982; Bardeen et al. 1983; Mukhanov et al. 1992). During the 
accelerated expansion phase, generation and amplification of quantum fiuctuations in scalar 
fields are unavoidable (Parker 1969; Birrell & Davies 1982). These fiuctuations become 
classical after crossing the event horizon. Later during the deceleration phase they re-enter 
the horizon, and seed the matter and the radiation fiuctuations observed in the universe. 

The majority of infiation models predict Gaussian, adiabatic, nearly scale-invariant 
primordial fiuctuations. These properties are generic predictions of infiationary models. The 
cosmic microwave background (CMB) radiation anisotropy is a promising tool for testing 
these properties, as the linearity of the CMB anisotropy preserves basic properties of the 
primordial fiuctuations. In companion papers, Spergel et al. (2003) find that adiabatic 
scale-invariant primordial fluctuations flt the WMAP CMB data as well as a host of other 
astronomical data sets including the galaxy and the Lyman-a power spectra; Komatsu et al. 
(2003) flnd that the WMAP CMB data is consistent with Gaussian primordial fluctuations. 
These results indicate that predictions of the most basic inflationary models are in good 
agreement with the data. 

While the inflation paradigm has been very successful, radically different inflationary 
models yield similar predictions for the properties of fluctuations: Gaussianity, adiabaticity, 
and near-scale-invariance. To break the degeneracy among the models, we need to measure 
the primordial fiuctuations precisely. Even a slight deviation from Gaussian, adiabatic, near- 
scale- invariant fiuctuations can place strong constraints on the models (Liddle & Lyth 2000). 
The CMB anisotropy arising from primordial gravitational waves can also be a powerful 
method for model testing. In this paper, we confront predictions of various inflationary 
models with the CMB data from the WMAP, CBI (Pearson et al. 2002), and ACBAR (Kuo 
et al. 2002) experiments, as well as the 2dFGRS (Percival et al. 2001) and Lyman-a power 
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spectra (Croft et al. 2002; Gnedin & Hamilton 2002). 

This paper is organized as follows. In § 2, we show that the WMAP detection of an 
anti-correlation between the temperature and the polarization fluctuations at I ~ 150 is 
the distinctive signature of adiabatic superhorizon fluctuations. We compare the data with 
specific predictions of inflationary models: single-fleld models in § 3, and double-field models 
in § 4. We examine the evidence for features in the inflaton potential in § 5. Finally, we 
summarize our results and draw conclusions in § 6. 

2. IMPLICATIONS OF WMAP "TE" DETECTION FOR THE 

INFLATIONARY PARADIGM 

A fundamental feature of inflationary models is a period of accelerated expansion in the 
very early universe. During this time, quantum fluctuations are highly amplified, and their 
wavelengths are stretched to outside the Hubble horizon. Thus, the generation of large-scale 
fluctuations is an inevitable feature of inflation. These fluctuations are coherent on what 
appear to be superhorizon scales at decoupling. Without accelerated expansion, the causal 
horizon at decoupling is ~ 2 degrees. Causality implies that the correlation length scale for 
fluctuations can be no larger than this scale. Thus, the detection of superhorizon fluctuations 
is a distinctive signature of this early epoch of acceleration. 

The COBE DMR detection of large scale fluctuations has been sometimes described as 
a detection of superhorizon scale fluctuations. While this is the most likely interpretation 
of the COBE results, it is not unique. There are several possible mechanisms for generating 
large-scale temperature fluctuations. For example, texture models predict a nearly scale- 
invariant spectrum of temperature fluctuations on large angular scales (Pen et al. 1994). 
The COBE detection sounded the death knell for these particular models not through its 
detection of fluctuations, but due to the low amplitude of the observed fluctuations. The 
detection of acoustic temperature fluctuations is also sometimes evoked as the definitive 
signature of superhorizon scale fluctuations (Hu & White 1997). String and defect models 
do not produce sharp acoustic peaks (Albrecht et al. 1996; Turok et al. 1998). However, 
the detection of acoustic peaks in the temperature angular power spectrum does not prove 
that the fluctuations are superhorizon, as causal sources acting purely through gravity can 
exactly mimic the observed peak pattern (Turok 1996a,b). The recent study of causal seed 
models by Durrer et al. (2002) shows that they can reproduce much of the observed peak 
structure and provide a plausible flt to the pre- WMAP CMB data. 

The large-angle (50 < I < 150) temperature-polarization anti-correlation detected by 
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WMAP (Kogut et al. 2003) is a distinctive signature of super horizon adiabatic fluctua- 
tions (Spergel & Zaldarriaga 1997). The reason for this conclusion is explained as follows. 
Throughout this section, we consider only scales larger than the sound horizon at the decou- 
pling epoch. Zaldarriaga & Harari (1995) show that, in the tight coupling approximation, 
the polarization signal arises from the gradient of the peculiar velocity of the photon fluid, 

ei, 

Ae ~ -0.17(1 - i2^)Ar]aeckei{r]dec), (1) 

where is the E-mode (parity-even) polarization fluctuation, rjdec is the conformal time 
at decoupling, Arjaec is the thickness of the surface of last scattering in conformal time, and 
fj, — cos{k ■ h). The velocity gradient generates a quadrupole temperature anisotropy pat- 
tern around electrons which, in turn, produces the £'-mode polarization. Note that while 
reionization violates the assumptions of tight coupling, the existence of clear acoustic os- 
cillations in the temperature-polarization (TE) and temperature-temperature (TT) angular 
power spectra imply that most (~ 85%) CMB photons detected by WMAP did indeed come 
from z = 1089 where the tight coupling approximation is valid. The velocity Gi is related to 
the photon density fluctuations, ©o, through the continuity equation, kQi = —3 ^Oq + 
where $ is Bardeen's curvature perturbation. The observable temperature fluctuations on 
large scales are approximately given by = ©o(%ec) + ^(%ec); where is the Newtonian 
potential, which equals — $ in the absence of anisotropic stress. Therefore, roughly speak- 
ing, the photon density fluctuations generate temperature fluctuations, while the velocity 
gradient generates polarization fluctuations. 

The tight coupling approximation implies that the baryon photon fluid is governed by 
a single second-order differential equation which yields a series of acoustic peaks (Peebles & 
Yu 1970; Hu & Sugiyama 1995): 

(©0 + ^) + ^ YT^(©o + ^)+ k'cKeo + ^)^k' (^cl^ - I) , (2) 

where the sound speed Cg is given by = [3(1 -|- i?)]~\ The large-scale solution to this 
equation is (Hu & Sugiyama 1995) 

60(77) + Hv) = [e>o(0) + $(0)] cos(A;c,r7) + fee, P drj' [^(ry') - ^{rj')] sm[kcs{v - V% (3) 







and the continuity equation gives the solution for the peculiar velocity, 

-^81(77) = [eo(0) + $(0)] sm{kcs7j) - kcs r dri' [$(V) - ^(V)] cos[kcs{ii - r/')]- (4) 

These solutions (equations (1), (3), and (4)) are valid regardless of the nature of the source 
of fluctuations. 
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In inflationary models, a period of accelerated expansion generates superhorizon adia- 

batic fluctuations, so that the first term in equation (3) and (4) is non-zero. Since ^ ~ 
— $ and 6o(0) + $(0) = |$(0) = |$(r7dec) on superhorizon scales, one obtains — 
-^^iVdec) cos{kCsr]dec), and A^^ ~ 0.17{l-fi'^)kCsAr]dec^{r]dec) sm{kcsr]dec) (see Hu & Sugiyama 
(1995) and Zaldarriaga & Harari (1995) for derivation). Therefore, the cross correlation is 
found to be 

(AtAe) ~ -0.03(1 - fi^){kcsAr]d,,)P^{k) sin(2A;c,%ec), (5) 

where P<i,{k) is the power spectrum of ^{ridec)- The observable correlation function is esti- 
mated as /c'^(ArAE). Clearly, there is an anti-correlation peak near kcsTjdec ~ 37r/4, which 
corresponds to I ~ 150: this is the distinctive signature of primordial adiabatic fluctua- 
tions. In other words, the anti-correlation appears on superhorizon scales at decoupling, 
because of the modulation between the density mode, cos(A;c577(iec), and the velocity mode, 
sm.{kcsr]dec)-i yielding sm.{2kcsr]dec)-i which has a peak on scales larger than the horizon size. 

Cosmic strings and textures are examples of active models. In these models, causal field 
dynamics continuously generate spatial variations in the energy density of a field. Magueijo 
et al. (1996) describe the general dynamics of active models. These models do not have the 
first term in equation (3) and (4), but the fluctuations are produced by the second term, 
the growth of $ and \E'. The same applies to primordial isocurvature fluctuations, where the 
non-adiabatic pressure causes $ and ^ to grow. While the problem is more complicated, 
these models give a positive correlation between temperature and polarization fluctuations 
on large scales. This positive correlation is predicted not just for texture (Seljak et al. 1997) 
and scaling seed models (Durrer et al. 2002), but is the generic signature of any causal models 
(Hu & White 1997)^^ that lack a period of accelerated expansion. 

Figure 1 shows the predictions of the TE large angle correlation predicted in typical 
primordial adiabatic, isocurvature, and causal scahng seed models compared with the WMAP 
data. The causal scaling seed model shown is a flat Family I model in the classification of 
Durrer et al. (2002) that provided a good fit to the pre- WMAP temperature data. 

The WMAP detection of a TE anti-correlation at / ~ 50 — 150, scales that correspond 
to superhorizon scales at the epoch of decoupling, rules out a broad class of active models. 
It implies the existence of superhorizon, adiabatic fiuctuations at decoupling. If these fiuc- 
tuations were generated dynamically rather than by setting special initial conditions then 
the TE detection requires that the universe had a period of accelerated expansion. In addi- 
tion to infiation, the pre-Big-Bang scenario (Gasperini & Veneziano 1993) and the Ekpyrotic 



Hu & White (1997) use an opposite sign convention for the TE cross power spectrum. 
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scenario (Khoury et al. 2001, 2002) predict the existence of superhorizon fluctuations. 

3. SINGLE FIELD INFLATION MODELS 

In this section we explore how predictions of specific models that implement inflation 
(see Lyth & Riotto (1999) for a survey) compare with current observations. 

3.1. Introduction 

The definition of "single-field infiation" encompasses the class of models in which the 
inflationary epoch is described by a single scalar field, the infiaton field. We also include a 
class of models called "hybrid" inflation models as single-fleld models. While hybrid inflation 
requires a second fleld to end inflation (Linde 1994), the second held does not contribute 
to the dynamics of inflation or the observed fluctuations. Thus, the predictions of hybrid 
inflation models can be studied in the context of single-fleld models. 

During inflation the potential energy of the inflaton fleld V dominates over the kinetic 
energy. The Priedmann equation then tells us that the expansion rate, H, is nearly constant 
in time: H = a/a M'^V/Sy/'^, where Mp, = (87rG')-V2 ^ rUpi/VSn = 2.4 x lO^^ GeV 
is the reduced Planck energy. The universe thus undergoes an accelerated expansion phase, 
expanding exponentially as a(t) oc exp(J Hdt) ~ exp{Ht). One usually uses the e-folds 
remaining at a given time, N{t), as a measure of how much the universe expands from t 
to the end of infiation, tend- N{t) = ln[a(tcnd)] — ln[a(t)] = H{t)dt. It is known that 
fiatness and homogeneity of the universe require A^(tstart) > 50, where tgtart is the time at 
the onset of inflation (i.e., the universe needs to be expanded to at least ~ 5 x 10^^ times 
larger by tend)- The accelerated expansion of this amount dilutes any initial inhomogeneity 
and spatial curvature until they become negligible in the observable universe today. 

3.2. Framework for data analysis 

3.2.1. Parameterizing the primordial power spectra 

The power spectrum of the CMB anisotropy is determined by the power spectra of 
the curvature and tensor perturbations. Most inflationary models predict scalar and tensor 
power spectra that approximately follow power laws: l\\{k) = A;V(27^2)(|7^k^ oc fc^--^ and 
Al{k) = 2A;^/(27r^)(|/i+kP + |^xk|^) oc k'^*. Here, TZ is the curvature perturbation in the 
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comoving gauge, and h+ and hx are the two polarization states of the primordial tensor 
perturbation. The spectral indices Ug and rit vary slowly with scale, or not at all. As 
spectral indices deviate more and more from scale invariance (i.e., = 1 and = 0), the 
power-law approximation usually becomes less and less accurate. Thus, in general, one must 
consider the scale dependent "running" of the spectral indices, dus/dhik and dnt/dlnk. We 
parameterize these power spectra by 

/ 7 \ ns{ko)—l+^{dna/dlnk)ln{k/ko) 

Alik) = Aliko)!-] , (6) 



Alik) = Al{ko)(^] , (7) 



^ \ nt{ko}+^ (dnt/dlnk) ln{k/ko} 



where A^(/co) is a normalization constant, and ko is some pivot wavenumber. The running, 
dn/dln k, is defined by the second derivative of the power spectrum, dn/ dink = ^A^j d In /c^, 
for both the scalar and the tensor modes, and is independent of k. This parameterization 
gives the definition of the spectral index, 

^ din A?, , , dus , f k\ 

for the scalar modes, and 

(i In A? , , , drit , f k\ 

for the tensor modes. In addition, we re-parameterize the tensor power spectrum amplitude, 
A^(A;o), by the "tensor/scalar ratio r", the relative amplitude of the tensor to scalar modes, 
given by^^ 

r = ^ (10) 

The ratio of the tensor quadrupole to the scalar quadrupole, r2, is often quoted when referring 
to the tensor/scalar ratio. The relation between r2 and the definition of the tensor/scalar 
ratio above is somewhat cosmology-dependent. For an SCDM universe with no reionization, 
it is: 

rs = 0.8625 r. (11) 

For comparison, for the maximum likelihood single field infiation model for the WMAPcxt-|-2dFGRS 
data sets presented in the table notes of Table 1 in §3.3, this relation is r2 = 0.6332 r. 



^^This definition of r agrees with the definition of T/S in the CAMB code (Lewis et al. 2000) and r in 
Leach et aL (2002). We have modified CMBFAST (Seljak & Zaldarriaga 1996) accordingly to match the same 
convention. 
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FoUowing notational conventions in Spergel et al. (2003), we use A{ko) for the scalar 
power spectrum amplitude, where A{ko) and A^{ko) are related through 

Aliko) = soon' -^A{ko) (12) 

\ "J / CMB 

~ 2.95 X 10-^^(A;o). (13) 

Here, Tcmb = 2.725 x 10^ (iiK). This relation is derived in Verde et al. (2003). One can 
use equations (6), (8), and (9) to evaluate A, Us, and rit at a different wavenumber from ko, 
respectively. Hence, 

/ 7 \ ne{ko) — l+^{dns/d\nk)ln{ki/ko) 

A{k,) = A{ko) [^J . (14) 

We have 6 observables {A, r, Ug, rit, dUg/dlnk, drit/dlnk), each of which can be com- 
pared to predictions of an inflationary model. 

The complementary approach (which we do not investigate in this work) is to parame- 
terize the primordial power spectrum in a model-independent way (see, for example, Wang 
et al. (1999)). These authors anticipated that WMAP has the potential ability to reveal 
deviations from scalc-invariance when combined with large scale structm^c data. Mukherjee 
& Wang (2003a,b) extend this approach and use it to put model-independent constraints on 
the primordial power spectrum using the pre- WMAP CMB data. 



3.2.2. Slow roll parameters 

In the context of slow roll inflationary models, only three "slow-roll parameters" , plus the 
amplitude of the potential, determine the six observables {A, r, Ug, fit, dus/dlnk, drit/dlnk). 
Thus, one can use the relations among the observables to either reduce the number of pa- 
rameters to four, or cross-check if the slow roll inflation paradigm is consistent with the data. 
The slow-roll parameters are defined by (Liddle & Lyth 1992, 1993): 



V'V- 



III 



^ ^pi(-^), (17) 



-9- 



where prime denotes derivatives with respect to the field 0. Here, quantifies "steepness" of 
the slope of the potential which is positive-definite, rjy quantifies "curvature" of the potential, 
and (which is not positive-definite, but is unfortunately often denoted in the literature 
because it is a second order parameter), quantifies the third derivative of the potential, or 
"jerk". All parameters must be smaller than one for infiation to occur. We denote these 
"potential slow roll" parameters with a subscript V to distinguish them from the "Hubble 
slow roll" parameters of Appendix A. Gratton et al. (2003) discuss the equivalent set of 
parameters for the Ekpyrotic scenario. 

Parameterization of slow roll models by ey, rjv, and avoids relying on specific models, 
and enables one to explore a large model space without assuming a specific model. Each 
inflation model predicts the slow-roll parameters, and hence the observables. A standard 
slow roll analysis gives observable quantities in terms of the slow roll parameters to flrst 
order as (see Liddle & Lyth (2000) for a review), 

r = l&ey, (19) 
Us -I = -Qey + 2r]y = -— + 27)y, (20) 

o 

nt = -2ey^~, (21) 

dv 3 2 

IQeyTiy - 2^el - 2^y = r7]y - -r^ - 2^y = -- [{ris - If - 44] - 2^^,(22) 



dink ^ 32 3 

drit 



dink 



AeyTiy - Se'y = - (n,-l) + ^ . (23) 



The tensor tilt in infiation is always red, Ut < 0. The equation Ut = — r/8 is known as the 
consistency relation for single-field infiation models (it weakens to an inequality for multi- 
field infiation models). We use the relation to reduce the number of parameters. While 
we have also carried out the analysis including rit as a parameter, and verified that there 
is a parameter space satisfying the consistency relation, including rit obviously weakens the 
constraints on the other observables. Given that we find r is consistent with zero (§ 3.3), 
the running tensor index drif/dlnk is poorly constrained with our data set; thus, we ignore 
it and constrain our models using the other four observables {A, r, Us, dus/dlnk) as free 
parameters. 
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3.3. Determining the power spectrum peirameters 

We use a Markov Chain Monte Carlo (MCMC) technique to explore the likelihood 
surface. Verde et al. (2003) describe our methodology. We use the WMAP TT (Hinshaw 
et al. 2003) and TE (Kogut et al. 2003) angular power spectra. To measure the shape of 
the spectrum (i.e., Ug and dng/dlnk) accurately, we want to probe the primordial power 
spectrum over as wide a range of scales as possible. Therefore, we also include the CBI 
(Pearson et al. 2002) and ACBAR (Kuo et al. 2002) CMB data, Lyman a forest data (Croft 
et al. 2002; Gnedin & Hamilton 2002), and the 2dFGRS large-scale structure data (Percival 
et al. 2001) in our hkelihood analysis. We refer to the combined WMAP+CBI+ ACBAR 
data as WMAPex.t. 

In total, the single field inflation model is described by an 8-parameter model: 4 parame- 
ters for characterizing a Friedmann-Robertson- Walker universe (baryonic density Qbh"^, mat- 
ter density Vtmh"^- Hubble constant in units of 100 kms~^Mpc^^ h, optical depth r), and 4 pa- 
rameters for the primordial power spectra {A, r, n^, dus/dlnk). When we add 2dFGRS data, 
we need two further large-scale structure parameters, [3 and dp, to marginalize over the shape 
and the amplitude of the 2dFGRS power spectrum (Verde et al. 2003). We run MCMC with 
these eight {WMAP only model) or ten ( WMAPext+2dFGRS, WMAPext+2dFGRS+Ly a 
models) parameters in order to get our constraints. 

The priors on the model are: a flat universe, a cosmological constant equation of state 
for the dark energy, and a restriction of r < 0.3. 



Table 1. Parameters For Primordial Power Spectra: Single Field Inflation Model 



Paraiiiotor 




iri\MPoxt+2(iFGRS" 


iri\MPoxt+2(iFGRS+L>-maii o;^ 




Usiko = 0.002 Mpc"^) 
r{ko = 0.002 Mpc"^) 
dus/dlnk 

A{ko = 0.002 Mpc-^) 


< 1.28/0.81/0.47^ 

-0 077+° °^° 

'-0.052 

71+010 

^-0.11 

0.024 ± 0.002 


1 1 C+O-12 

< 1.14/0.53/0.37'' 
-0 075+°-°^^ 
0.73 ± 0.09 
0.023 ±0.001 


1.13 ±0.08 
< 0.90/0.43/0.29*^ 
-0 n55+°-°28 

n yr+O.OS 
^- ' ^-0.09 

0.024 ±0.001 






0.127 ±0.017 


0.134 ±0.006 


0.134 ±0.006 




h 


0.78 ±0.07 


' J-0.04 


0.75 ±0.03 




T 


0.22 ±0.06 


0.20 ±0.06 


0.18 ±0.06 


1 


(78 


o.82iS:i^ 


0.85 ±0.05 


0.85 ± 0.05 


I— » 
I— ' 



I 



'^The quoted values are the mean and the 68% probability level of the 1-d marginalized likelihood. For 
both WMAPext±2dFGRS and WMAPext±2dFGRS±Lyman a data sets, the 10-d maximum hkelihood 
point in the Markov Chain (1.5 x 10^ steps) for this model is [n^/i^ = 0.024, fi^/i^ ^ o.l32, h = 0.77, 
ri(A;o.oo2) = 1-15, r(A;o.oo2) = 0.42, dus/dlnk = -0.052, ^(W) = 0.75, r = 0.21, erg = 0.87]. Here, 
^0.002 is /co = 0.002 Mpc~^. The maximum hkelihood model in the MCMC using WMAP data alone 
is [Qfe/i^ ^ 0.023, Qrnh'^ = 0.122, h = 0.79, n(A;o.oo2) = 1-27, r(A;o.oo2) = 0.56, dus/dlnk = -0.10, 
^(^0.002) = 0.74, r = 0.29]. Great care must be taken in interpreting this point. It is given here for 
completeness only, and we do not recommend it for use in any analysis. There is a long, fiat degeneracy 
between n and r, as described in §3 Spergel et al. (2003), and this point happened to lie at the very 
blue edge of this degeneracy right at the edge of our upper limit prior on r. This Markov chain had 
extra freedom because we are adding three parameters over the model discussed in Spergel et al. (2003), 
thereby introducing significant new degeneracies (see Figure 3). 

^The 95% upper limits for the tensor-scalar ratio are quoted for various priors in the following order: 
[no prior on dUs/dlnk or Ug] / [dng/dlnk— 0] / [ug < 1]. The priors were applied to the output of the 
MCMC. 
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Table 1 shows results of our analysis for the WMAP, WMAPex.t+ 2dFGRS and WMAPext+ 
2dFGRS + Lyman a data sets. We evaluate Ug, A, and r in the fit at ko = 0.002 Mpc~^ 
Thus, this table and the figures to follow report the results for A and Ug at ko = 0.002 Mpc^^ 
. Note that Spergel et al. (2003) report these quantities evaluated at ko = 0.05 Mpc^^ 
(using equations (14) and (8)). There are 3.2 e-folds between ko = 0.002 Mpc~^ and 
ko = 0.05 Mpc-^ 

We did not find any tensor modes. Table 1 shows 95% upper limits for the tensor- 
scalar ratio r at k = 0.002 Mpc~^, for various combinations of the data sets. As we will see 
later, there arc strong degeneracies present between the parameters n.^, r and dus/dlnk. For 
example, one can add power at low multipoles by increasing r and then remove it with a bluer 
Us while keeping the low / amplitude constant. Thus, one can obtain stronger constraints 
on r by assuming different priors on Ug and dus/dlnk. In the table, we list the 95% CL 
constraints on r that would be obtained if (1) there were no priors on Ug or dus/dlnk, (2) if 
one only considers models with no running of the scalar spectral index, and (3) if only models 
with red spectral indices are considered (non-hybrid-inflation models predict red indices in 
general) . 

The no-prior r limit r < 0.9, along with the 2-a upper limit on the amplitude A{k = 
0.002 Mpc"^) < 0.75 + 0.08 x 2, implies that the energy scale of inflation F V4 < 3.3 x 10^^ 
GeV at the 95% confidence level. 

Note that in the case of the WMAP-only Markov chain, the degeneracy between Ug, 
r and dug/dlnk is cut off by the prior r < 0.3 (r is denegerate with n^). Thus, a better 
upper limit on r will significantly tighten the constraints on this model from the WMAP 
data alone. 

All cosmological parameters are consistent with the best-fit running model of Spergel 
et al. (2003), which was obtained for a ACDM model with no tensors and a running spectral 
index. Adding the extra parameter r does not improve the fit. 

Our constraint on Ug shows that the scalar power spectrum is nearly scale- invariant. One 
implication of this result is that fluctuations were generated during accelerated expansion 

in nearly de-Sitter space (Mukhanov & Chibisov 1981; Hawking 1982; Guth & Pi 1982; 
Starobinsky 1982; Bardcen et al. 1983; Mukhanov et al. 1992), where the equation of state 
of the scalar field is w ~ — 1. Recently, Gratton et al. (2003) have shown that there is only 
one other possibility for robustly obtaining adiabatic fluctuations with nearly scale-invariant 
spectra: w 3> 1. The Ekpyrotic/ Cyclic scenarios correspond to this case. Note, however, that 
predictions for the primordial perturbation spectrum resulting from the Ekpyrotic scenario 
are controversial (see, for example, Tsujikawa et al. (2002)). 
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We find a marginal 2a preference for a running spectral index in all three data sets; 
dn^/dlnfe = -0.055l2;|]29 ( WMAPext+2dFGRS+Lyman a data set). This same preference 
was seen in the analysis without tensors carried out in Spergel et al. (2003). 

Figure 2 shows our constraint on function of k for the WMAF, WMAF^-yiX^ 

2dFGRS and WMAF^sX^ 2dFGRS + Lyman a. data sets. At each wavenumber A;, we use 
equation (8) to convert ns(A;o = 0.002 Mpc~^) to risili) at each wavenumber. Then, we 
evaluate the mean (solid line), 68% interval (shaded area), and 95% interval (dashed lines) 
from the MCMCs. This shows a hint that the spectral index is running from blue [rig > 1) 
on large scales to red (n^ < 1) on small scales. In our MCMCs, for the WMAP data set 
alone, 91% of models explored by the chain have a scalar spectral index running from blue 
at A; = 0.0007 Mpc"^ (/ ~ 10) to red at A; = 2 Mpc'^ For the WMAPext+2dFGRS data 
set, 95% of models go from a blue index at large scales to a red index at small scales, and 
when Lyman a forest data is added, the fraction running from blue to red becomes 96%. 

One-loop correction and renormalization usually predict running mass and/or running 
coupling constant, giving some dus/dlnk. Detection of it implies interesting quantum phe- 
nomena during inflation (see Lyth & Riotto (1999) for a review). For the running of the 
scalar spectral index (equation 22), 

J?^ = -2«.-H[K-I)=-4,a. (24) 

Since the data requires ~ 1 (see Table 1), {ug — 1)^ < 0.01. It is especially small when 
— 1 ~ 2r]v, (see Case A and Case D in § 3.4.2). Therefore, if dus/dlnk is large enough to 
detect, dus/dlnk > 10"^, then dug/dlnk must be dominated by 2^y, a product of the first 
and the third derivatives of the potential (equation (17)). The hint of dUs/dlnk in our data 
can be interpreted as — —-^dns/dlnk = 0.028 ± 0.015. However, obtaining the running 
from blue to red, which is suggested by the data, may require fine-tuned properties in the 
shape of the potential. More data are required to determine whether the hints of a running 
index are real. 



3.4. Single field models confront the data 

3.4- i- Testing a specific inflation model: \(f)^ 

As a prelude to showing constraints on broad classes of inflationary models, we first 
illustrate the power of the data using the example of the minimally-coupled V — A0^/4 
model, which is often used as an introduction to infiationary models (Linde 1990). We show 
that this textbook example is unlikely. 
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The Priedmann and continuity equations for a homogeneous scalar field lead to the 
slow-roll parameters, which one can use in conjunction with the equations of § 3.2.2 in order 
to obtain predictions for the observables. For the potential 1^(0) = A0^/4, one obtains the 
potential slow roll parameters as: 

Ml, M2 M^i 

The number of e-foldings remaining till the end of infiation is defined by 
N 



where ty{(j)f,nd) = 1 defines the end of inflation. Assuming (p^nd ^ (p, taking the horizon exit 
scale as ~ V SNMpi and N — 50, one obtains Ug — 0.94 and r = 0.32 using equations (19) 
and (20). As dns/dlnk is negligible for this model, we use dus/dlnk — 0. 

We maximize the likelihood for this model by running a simulated annealing code. 
We fit to WMAPext -|-2dFGRS data, varying the following parameters: Q.f,h'^, ^m.h'^: h, r, 
A^^, (3, and Up, while keeping Ug, dus/dlnk, and r fixed at the X(f)^ values. The maximum 
likehhood model obtained has [Q^/i^ = 0.022, Q^/i^ ^ o.l35, r = 0.07, A = 0.67, h = 
0.69, CTg = 0.76]. This best-fit model is compared in Table 2 to the corresponding model 
with the full set of single field inflationary parameters. The model is displaced from 
the maximum likelihood generic single field model by ^xlff = 16 [Axlff{WMAP) = 14, 
AXe//(CBI+ACBAR+2dFGRS) = 2], where xljf = -21n£ and C is the likelihood (see 
Verde ct al. (2003)). Since the relative likelihood between the models is cxp(— 8), and the 
number of degrees of freedom is approximately three, A^"^ is disfavored at more than 3a. 
The table shows that adding external data sets does not make a significant difference to the 
Axgjj between the models, and the constraint is primarily coming from WMAP data. 

This result holds only for Einstein gravity. When a non-minimal coupling of the form 
i(f)^R (.^ = 1/6 is the conformal coupling) is added to the Lagrangian, the coupling changes 
the dynamics of 0. This model predicts only a tiny amount of tensor modes (Komatsu & 
Futamase 1999; Hwang & Noh 1998) in agreement with the data. 



^^While A is an inflationary parameter, it is directly related to the self-coupling A which we do not know; 
thus, we treat it as a parameter. 
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Table 2. Goodness-of-fit Comparison for Model 



Model 


Xlff (WMAP) 


xlff (ext+2dFGRS) 


Total xlff/ 1' (WMAPext+2dFGRS) 


Best-fit inflation 


1428 


36 


1464/1379 


A0*^ model 


1442 


38 


1480/1382 
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One can perform a similar analysis on any given inflationary model to see what con- 
straints the data put on it. Rather than attempt this Herculean task, in the following section 
we simply use our constraints on n^, dns/dlnk, and r and the predictions of various classes 
of single field inflationary models for these parameters in order to put broad constraints on 
them. 



3.4-2. Testing a broad class of inflation models 

Naively, the parameter space in observables spanned by the slow roll parameters appears 
to be large. We shall show below that "viable" slow roll inflation models (i.e. those that can 
sustain inflation for a sufficient number of e-folds to solve cosmological problems) actually 
occupy significantly smaller regions in the parameter space. 

Hoffman & Turner (2001); Kinney (2002a); Easthcr & Kinney (2002); Hansen & Kunz 
(2002); Caprini et al. (2003) have investigated generic predictions of slow roll inflation models 
by using a set of inflationary flow equations (see Appendix A for a detailed description and 
deflnition of conventions). In particular Kinney (2002a) and Easther & Kinney (2002) use 
Monte Carlo simulations to extend the slow roll approximations to fifth order. These authors 
find "attractors" corresponding to fixed points (where all derivatives of the flow parameters 
vanish); models cluster strongly near the power-law inflation predictions, r — 8{1 — Us) (see 
§ 3.4.4), and on the zero tensor modes, r — 0. 

Following the method of Kinney (2002a) and Easther & Kinney (2002), we compute 
a million realizations of the inflationary flow equations numerically, truncating the flow 
equation hierarchy at eighth order and evaluating the observables to second order in slow 
roll using equations (A15)-(A17). We marginahze over the ambiguity of converting between 
(j) and k, introduced by the details of reheating and the energy density during inflation by 
adopting the Monte Carlo approach of the above authors. The observable quantities of 
a given realization of the flow equations are evaluated at a specific value of e-folding, A^. 
However, observable quantities are measured at a specific value of k. Therefore, we need 
to relate N to k. This requires detailed modeling of reheating, which carries an inherent 
uncertainty. We attempt to marginalize over this by randomly drawing A^ values from a 
uniform distribution A" = [40, 70]. 

Figure 3 shows part of the parameter space of viable slow roll infiation models, with the 
WMAP 95% confidence region shown in blue. Each point on these panels is a different Monte 
Carlo realization of the flow equations, and corresponds to a viable slow roll model. Not all 
points that are viable slow roll models correspond to speciflc physical models constructed in 
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the literature. Most of the models cluster near the attractors, sparsely populating the rest 
of the large parameter space allowed by the slow roll classification. It must be emphasized 
that these scatter plots should not be interpreted in a statistical sense since we do not know 
how the initial conditions for the universe are selected. Even if a given realization of the 
flow equations does not sit on the attractor, this does not mean that it is not favored. Each 
point on this plot carries equal weight, and each is a viable model of inflation. Notice that 
the WMAP data do not lie particularly close to the r — S{1 — Ug) "attractor" solution, at 
the 2-0" level, but is quite consistent with the r — attractor. 

One may categorize slow roll models into several classes depending upon where the 
predictions lie on the parameter space spanned by n^, dus/dlnk, and r (Dodelson et al. 
1997; Kinney 1998; Hannestad et al. 2001). Each class should correspond to specific physical 
models of infiation. Hereafter, we drop the subscript V unless there is an ambiguity — 
it should otherwise be implicitly assumed that we are referring to the standard slow roll 
parameters. We categorize the models on the basis of the curvature of the potential rj, as 
it is the only parameter that enters into the relation between Ug and r (equation (20)), and 
between Ug and dns/d\nk + 2^ (equation (22)). Thus, rj is the most important parameter 
for classifying the observational predictions of the slow roll models. The classes are defined 
by 

(A) negative curvature models, rj < 0, 

(B) small positive (or zero) curvature models, < 77 < 2e, 

(C) intermediate positive curvature models, 2e < 77 < 3e, and 

(D) large positive curvature models, rj > 3e. 

Each class occupies a certain region in the parameter space. Using r] = [ug — 1)q;/[2(q; — 3)], 
where 77 = cte, one finds 

(A) ris < 1, < r < 1(1 - Ug), -|(1 - Ugf < dng/d\nk + 2^ < 0, 

(B) Ug < 1, 1(1 -ng)<r< 8(1 - n,), -|(1 - Ugf < dng/dlnk + 2^ < 2(1 - Ugf, 

(C) 7^s < 1, r > 8(1 - Ug), dng/d\nk + 2^ > 2(1 - Ugf, and 

(D) > 1, r > 0, dng/dlnk + 2{ > 0. 

To first order in slow roll, the subspace {ug, r) is uniquely divided into the four classes, and 
the whole space spanned by these parameters is defined by this classification. The division 
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of the other subspace (n^, dus/lnk) is less unique, and dus/dlnk < —2^ — |(1 — Ug)^ is 
not covered by this classification. To higher order in slow roll, these boundaries only hold 
approximately - for instance. Case C can have a slightly blue scalar index, and Case D can 
have a slightly red one. 

We summarize basic predictions of the above model classes to first order in slow roll 
using the relation between r and Us (equation (20)) rewritten as 

r = ^(l-n,) + y77. (27) 

This imphes: 

(A) negative curvature models predict 77 < and 1 — > 0; the second term nearly cancels 
the first to give r too small to detect, 

(B) small positive curvature models predict 1 — > and ?] > 0; a large r is produced, 

(C) intermediate positive curvature models predict 1 — Ug > and 77 > 0; a large r is 
produced, and 

(D) large positive curvature models predict 1 — < and 77 > 0; the first term nearly 
cancels the second to give r too small to detect. 

The cancellation of the terms in Case A and Case D implies — 1 ~ 277: the steepness 
of the potential in Case A and Case D is insignificant compared to the curvature, e ^ |77|. 
On the other hand, in Case B and Case C the steepness is larger than or comparable to 
the curvature, by definition; thus, non-detection of r can exclude many models in Case B 
and Case C. As we have shown in § 3.4.1, a minimally-coupled A0'^ model, which falls into 
Case B, is excluded at high significance, largely because of our non-detection of r (see also 
§ 3.4.4). 

For an overview. Figure 4 shows the Monte Carlo fiow equation realizations correspond- 
ing to the model classes A-D above on the {ug, r), {ug, dUg/dln k), and (r, dUg/dln k) planes, 
for the WMAP, WMAPext+2dFGRS and WMAPext+2dFGRS+Lyman a data sets. 

In Table 3, we show the ranges taken by the obscrvablcs n^, r and dng/ In A; in the Monte 
Carlo realizations that remain after throwing out all the points which are outside at least one 
of the joint-95% confidence levels. These points have been separated into the model classes 
A-D via their 77^.. These constraints were calculated as follows. First, we find the Monte 
Carlo realizations of the fiow equations from each model class that fall inside all the joint- 
95% confidence levels for a given data set, separately for the WMAP, WMAPext-F2dFGRS 



-19- 



and WMAPext+2dFGRS+Lyman a data sets (i.e. the models shown on Figure 4). Then 
wc find for each model class the maximum and minimum values predicted for each of the 
observables within these subsets. These constraints mean that only those models (within 
each class) predicting values for the observables that lie outside these limits arc excluded by 
these data sets at 95% CL. Note that the best-fit model within this parameter space has a 
xlff/^ = 1464/1379. Here, recall again that the observables were evaluated to second order 
in slow roll in these calculations. This is the reason that the Class C range in Ug goes slightly 
blue and the Class D range in Ug goes shghtly red; the divisions of the rjy classification are 
only exact to first order in slow roll. 

In the following subsections we will discuss in more detail the constraints on specific 
physical models that fall into the classes A-D. For a given class, we will plot only the flow 
equation realizations falling into that category that are consistent with the 95% confidence 
regions of all the planes {us, r), {ug, dus/dlnk) and (r, dng/dlnk). 



Table 3. Properties of Inflationary Models Present Within the Joint-95% Confidence 



Region^ 


Model 


X T 7"Tl /r A "T"^ 

WMAP 


WMAPext+2dFGRS 


WMAPext+2dFGRS+Lyman a 


A 


(A X in"^l^< r < n 14 


(2 X 10^^1^< r < 19 


(A X 10~^1'^< r < Ifi 






93 < T7 < 1 00 


94 < T7 < 1 00 

\J . Urt _^ 1 tg ± .\J\J 




-0.02 < dus/dlnk < 0.02 


-0.04 < dus/dlnk < 0.02 


-0.02 < dus/dlnk < 0.004 


B 


(7 X 10-=^)^< r < 0.35 


(7 X 10-^)^< r < 0.32 


(7 X 10-=^)^< r < 0.26 




0.94 < < 1.01 


0.93 < < 1.01 


0.94 <ns< 1.01 




-0.02 < dnjdlnk < 0.02 


-0.04 < dus/dlnk < 0.02 


-0.02 < dus/dlnk < 0.01 


C 


(0.003)^< r < 0.59 


(0.003)^< r < 0.52 


(0.03)^< r < 0.46 




0.95 < < 1.02 


0.96 < n, < 1.02 


0.97 <ns< 1.02 




-0.04 < dUs/dlnk < 0.01 


-0.04 < dnjdlnk < 0.01 


-0.04 < dris/dXnk < 0.001 


D 


0.0 < r < 1.10 


0.0 < r < 0.89 


(8 X 10-5)^< r < 0.89 




0.99 < < 1.28 


1.00 < < 1.28 


1.00 < < 1.28 




-0.09 < dus/dlnk < 0.03 


-0.09 < dus/dluk < 0.01 


-0.09 < dus/dlnk < -0.001 



^The ranges taken by the predicted observables of slow roll models (to second order in slow roll) 
within the joint 95% CLs from the specified data sets. The model classes are: Case A (77 < 0), Case 
B (0 < 77 < 2e), Case C (2e < 77 < 3e), Case D (77 > 3e). 

•^The lower value of r does not represent a detection, but rather the minimal level of tensors predicted 
by any point in the Monte Carlo that falls within in this class and is consistent with the data. We 
include the lower limit to help set goals for future CMB polarization missions. 
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Note that very few models predict a "bad power law", or \dns/dlnk\ > 0.05. 

3.4-3. Case A: negative curvature models rj < 

The top row of Figure 5 shows the Monte Carlo points belonging to Case A which are 
consistent with all the joint-95% confidence regions of the observables shown in the figure, 
for the WMAPext+2dFGRS+Lyman a data set. 

The negative rj models often arise from a potential of spontaneous symmetry breaking 
(e.g., new inflation - Albrecht & Steinhardt (1982); Linde (1982)). 

Wc consider negative-curvature potentials in the form oi V = A"^[l — {(f) /fiY] where 
p > 2. Wc require (f) < fi for the form of the potential to be valid, and A determines the 
energy scale of inflation, or the energy stored in a false vacuum. One flnds that this model 
always gives a red tilt < 1 to first order in slow roll, as — 1 = — 6e — 2 |?7| < 0. 

For p = 2, the number of e-folds at (p before the end of inflation is given by ~ 
(/x^/2Mpj) ln(/x/0), where we have approximated 0end — yW- By using the same approximation, 
one finds - 1 ~ -4(Mpi//i)2, and r ~ 32(0^ M^jZ/x^) ~ 8(1 - nje-^^^""^). In this class 
of models, ris cannot be very close to 1 without /i becoming larger than Wpi. For example. 
Us — 0.96 imphes // ~ lOMpi ~ 2mpi. For this class of models, r has a peak value of r ~ 0.06 
at Us — 0.98 (assuming A^ = 50). Even this peak value is too small for WMAP to detect. We 
sec from Table 3 that this model is consistent with the current data, but requires /i > rupi 
to be valid. 

For p > 3, ris - 1 ~ -{2/N){p - l)/(p - 2) or 0.92 < < 0.96 for N ^ 50 regardless 
of a value of and r ~ 4p^(Mpi///)^(0///)^(^"^) is negligible as <^ These models he in 
the joint 2-a contour. 

The negative 77 model also arises from the potential in the form ofV — A^[l-|-Q;ln(0///)], 
a one-loop correction in a spontaneously broken supersymmetric theory (Dvali et al. 1994). 
Here the coupling constant a should be smaller than of order 1. In this model (f) rolls down 
towards the origin. One finds — 1 = — (1 + ^a)/N which implies 0.95 < Ug < 0.98 
for 1 > a > (this formula is not valid when a = or = yu). Since r = 8a /N = 
8a{l + |q;)"-'^(1 — Ug) = 0.016(«/0.1), the tensor mode is too small for WMAP to detect, 
unless the coupling a takes its maximal value, a ~ 1. This type of model is consistent with 
the data. 



-22- 



3.4-4- Case B: small positive curvature models < rj <2e 

The second row of Figure 5 shows the Monte Carlo points belonging to Case B which are 
consistent with all the joint-95% confidence regions of the observables shown in the figure. 

The "small" positive r) models correspond to monomial potentials for < rj < 2e 
and exponential potentials for t] = 2e. The monomial potentials take the form of V = 
A'^(0/yu)P where p >2, and the exponential potentials V = A^exp(0/yu). The zero rj model 
is V = A'^{(j)/n). To first order in slow roll, the scalar spectral index is always red, as 
Us — 1 — — 6e + 2r] < — 4e < 0. The zero rj model marks a border between the negative r] 
models and the positive rj models, giving r = |(1 — ris). 

The monomial potentials often appear in chaotic infiation models (Lindc 1983), which 
require that (p be initially displaced from the origin by a large amount, ~ mpi, in order to 
avoid fine-tuned initial values for 0. The monomial potentials can have a period of inflation 
at > THpi, and inflation ends when rolls down to near the origin. For p — 2, inflation 
is driven by the mass term, which gives = 2y/NMpi, Ug — 1 — 2/N — 0.96, r = = 
4(1 - Us) = 0.16, and dn^/dlnk = -2/N^ = -(1 - 71^)^/2 = -0.8 x 10"=*. For p = 4, 
inflation is driven by the self-coupling, which gives (j) = 2\/2NMp\, Ug = 1 — 3/N = 0.94, 
r = 16/A^ = f (1-n,) = 0.32, and dris/dlnk = -3/N^ = -(1 -n,)V3 = -1.2 x lO^^. The 
most striking feature of the small positive rj models is that the gravitational wave amplitude 
can be large, r > 0.16. Our data suggest that, for monomial potentials to lie within the 
joint 95% contour, r < 0.26 (Table 3). A model is excluded at ~ 3-a (§ 3.4.1), and any 
monomial potentials with p > 4 are also excluded at high signifcance. Models with p — 2 
(mass term inflation) are consistent with the data. 

The exponential potentials appear in the Brans-Dicke theory of gravity (Brans & Dicke 
1961; Dicke 1962) conformally transformed to the Einstein frame (the extended inflation 
models) (La & Steinhardt 1989). One flnds Ug — 1 — (/x/Mpi)^, r = 8(1—77-5), and dng/dlnk — 
0. Thus, the exponential potentials predict an exact power-law spectrum and signiflcant 
gravitational waves for signiflcantly tilted spectra. Since // = NM^^/{(f) — (f)end): rig — 1 — 
[NMpi/{(j) - (pendiV- The 95% range for rig in Table 3 implies that - 0end > 4A^Mpi ~ 
200Mpi ~ 40mpi. 

The exponential potentials mark a border between the small positive 77 models and the 
positive intermediate 77 models described below. 
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3.4-5. Case D: large positive curvature models rj > 3e 

Before describing Case C, it is useful to describe Case D first. The fourth row of Figure 
5 shows the Monte Carlo points belonging to Case D which are consistent with all the 
joint-95% confidence regions of the observables shown in the figure. 

The "large" positive curvature models correspond to hybrid infiation models (Linde 
1994), which have recently attracted much attention as an /^-invariant supersymmetric the- 
ory naturally realizes hybrid inflation (Copeland et al. 1994; Dvali et al. 1994). While it is 
pointed out that supergravity effects add too large an effective mass to the inflaton field to 
maintain inflation, the minimal Kahler supergravity does not have such a large mass problem 
(Copeland et al. 1994; Linde & Riotto 1997). The distinctive feature of this class of models 
with ?7 > 3e is that the spectrum has a blue tilt, — 1 = —66 + 27] > 0, to flrst order in slow 
roll. 

A typical potential is a monomial potential plus a constant term, V — A'^[l + (0///)^], 
which enables inflation to occur for a small value of 0, < mpi. At flrst sight, inflation never 
ends for this potential, as the constant term sustains the exponential expansion forever. 
Hybrid inflation models postulate a second fleld a which couples to 0. When rolls slowly 
on the potential, cr stays at the origin and has no effect on the dynamics. For a small value 
of inflation is dominated by a false vacuum term, V{(j), o" = 0) ~ A"^. When rolls down to 
some critical value, a starts moving toward a true vacuum state, V{(j), a) = 0, and inflation 
ends. A numerical calculation (Linde 1994) suggests that the potential is described by 
only until reaches a critical value. When reaches the critical value, inflation suddenly 
ends and a need not be considered. Thus, we include the hybrid models in our discussion of 
single-fleld models. 

For p = 2, one flnds that ~ |(/i./Mpi)^ ln(0/0end) — 50, which, in turn, implies 
fjL ~ lOMpi ~ 2mpi for ln(0/0end) ~ 1- The spectral slope is estimated as ~ 1 + 
4(Mpi///)2 ~ 1.04, and the tensor/scalar ratio, r ~ 32(0///)2(Mpi///)2 = 8((j)/ /^^(ns - 1), 
is negligible as inflation occurs at -C /x. The running is also negligible, as dus/dhik ~ 
64(0//x)2(Mpi//x)4 = 4(0//x)2(n, - If < 10-2. This type of model lies within the joint 95% 
contours. 

One-loop correction in a softly broken supersymmetric theory induces a logarithmically 
running mass, = A^ {1 -|- (0///)^ [1 + Q;ln(0/Q)]}, where a is a coupling constant and Q 
is a renormalization point. Since ris is practically determined by V", this potential gives 
rise to a logarithmic running of rig (Lyth & Riotto 1999). These models would lie in the 
region occupied by the Monte Carlo points that have a large, negative dus/dlnk. This type 
of model is consistent with the data. 
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3.4-6. Case C: intermediate positive curvature models 2e < rj < 3e 

The third row of Figure 5 shows the Monte Carlo points belonging to Case C which are 
consistent with all the joint-95% confidence regions of the observables shown in the figure. 

The "intermediate" positive curvature models are defined, to first order in slow roll, as 
having a red tilt, — 1 = — 6e + 2?] < 0, or the exactly scale-invariant spectrum, n., — 1 = 0, 
while not being described by monomial or exponential potentials. These conditions lead to 
a parameter space where 2e < rj < 3e. Here we discuss only examples of physical models 
that do not solely live in Case C, but briefly pass through it as they transition from Case D 
to Case B or Case A. 

The transition from Case D to Case B may correspond to a special case of hybrid 
inflation models described in the previous subsection (Case D), V = A''[l + {(p/n)^]. When 
(f) ^ fjL, the potential becomes Case B potential, V — > A^(0//i)^, and the spectrum is red, 
rig < 1. When (f) <^ /i, the potential drives hybrid inflation, and the spectrum is blue, rig > 1. 
On the other hand, when ~ the potential takes a parameter space somewhere between 
Case B and Case D, which corresponds to Case C. One may argue that this model requires 
fine-tuned properties in that we just transition from one regime to the other. However, the 
Case C regime has an interesting property: the spectral index Ug runs from red on large 
scales to blue on small scales, as undergoes the transition from Case B to Case D. This 
example has the wrong sign for the running of the index compared to the data at the ~ 2-a 
level. 

Linde & Riotto (1997) is one example of a transition from Case D to Case A. They 
consider a supergravity-motivatcd hybrid potential with a one-loop correction, which can be 
approximated during infiation as 

y ~ A^ [1 + a ln(0/g) + A(0//i)^] . (28) 

Suppose that the one-loop correction is negligible in some early time, i.e., ~ Q. The 
spectrum is blue. (The third term is practically unimportant, as infiation is driven by the 
first term at this stage.) If the loop correction becomes important after several e-folds, then 
rig changes from blue to red, as the loop correction gives a red tilt as we saw in § 3.4.3. This 
example is consistent with the data. The transition (from Case D to Case A) is possible 
only when a and Q conspire to balance the first term and the second term right at the scale 
accessible to our observations. 
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4. MULTIPLE FIELD INFLATION MODELS 
4.1. Framework 

In general, a candidate fundamental theory of particle physics such as a supersymmetric 
theory requires not only one, but many other scalar fields. It is thus naturally expected that 
during inflation there may exist more than one scalar field that contributes to the dynamics 
of inflation. 

In most single-field inflation models, the fluctuations produced have an almost scale- 
invariant, Gaussian, purely adiabatic power spectrum whose amplitude is characterized by 
the comoving curvature perturbation, TZ, which remains constant on superhorizon scales. 
They also predict tensor perturbations with the consistency condition in equation (21). 

With the addition of multiple fields, the space of possible predictions widens consider- 
ably. The most distinctive feature is the generation of entropy, or isocurvature, perturbations 
between one field and the other. The entropy perturbation, 5, can violate the conservation of 
iZ on superhorizon scales, providing a source for the late-time evolution of IZ which weakens 
the single field consistency condition into an upper bound on the tensor/scalar ratio (Po- 
larski & Starobinsky 1995; Sasaki & Stewart 1996; Garcia-Bellido & Wands 1996). Limits 
on the possible level of the entropy perturbation thus discriminate between the multiple field 
models and the single field models. In this section, we consider the minimal extension to 
single-field inflation - a model consisting of two minimally-coupled scalar fields. 

4.2. Correlated Adiabatic/Isocurvature Fluctuations from Double-Field 

Inflation 

The WMAP data confirm that pure isocurvature fluctuations do not dominate the 
observed CMB anisotropy. They predict large-scale temperature anisotropics that are too 
large with respect to the measured density fluctuations, and have the wrong peak/trough 
positions in the temperature and polarization power spectra (Hu & White 1996; Page et al. 
2003). The WMAP observations hmit but do not preclude the possibility of correlated 
mixtures of isocurvature and adiabatic perturbations, which is a generic prediction of two- 
fleld inflation models. Both isocurvature and adiabatic perturbations receive signiflcant 
contributions from at least one of the scalar flclds to produce the correlation (Langlois 1999; 
Pierpaoh et al. 1999; Langlois & Riazuelo 2000; Gordon et al. 2001; Bartolo et al. 2001, 2002; 
Amendola et al. 2002; Wands et al. 2002). We focus on these mixed models in this section. 

Let iZrad and Srad bc the curvature and entropy perturbations deep in the radiation era. 
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respectively. At large scales, the temperature anisotropy is given by (Langlois 1999): 

AT 



T 



— {^rad — '^Sradj , (29) 





in addition to the integrated Sachs-Wolfe effect. The entropy perturbation, Srad = Spcdm/ Pcdm- 
(3/4)5p^/p^, remains constant on large scales until re-entry into the horizon. If TZrad and 
Srad have the same sign (correlated), then the large scale temperature anisotropy is reduced. 
If they have opposite signs (anti-correlated), then the temperature anisotropy is increased. 
Spergel et al. (2003) find that there is an apparent lack of power at the very largest scales 
in the WMAP data; thus, one of the motivations of this study is to see whether a correlated 
Srad can provide a better fit to the WMAP \ow-l data than a purely adiabatic model. 

The evolution of the curvature/entropy perturbations from horizon-crossing to the 
radiation-dominated era can be parameterized by a transfer matrix (Amendola et al. 2002), 

(30) 

k=aH 

Here, Trr = 1 and Tsr = because of the physical requirement that TZ is conserved for 
purely adiabatic perturbations, and that TZ cannot source S. All the quantities in equation 
(30) are weakly scale-dependent, and may be parameterized by power-laws. Hence, we write 
this equation as 

^^rad ~~ ^ CL^ ~\~ A^k '^ftg, (^"^) 

Srad = Bk'^'ds, (32) 

where a-nd Og are independent Gaussian random variables with unit variance, {drdg) — Srs- 
The cross-correlation spectrum is given by A.\g{k) = {k^ / 27:'^) (TZradSrad) — ^s-B/c"2+"3_ Qne 
may define the correlation coefficient using an angle A as 

{TZradSrad) sign{B) A sk""^ 

COS A = — = — , (33) 

C^ldV^'i^ldV^' ./W^T^k^ 

where —1 < cos A < 1. Thus, in general, six parameters (A,., A,,, cos A, ni, n2, n^) are 
needed to characterize the double-inflation model with correlated adiabatic/isocurvature 
perturbations, while cos A is scale-dependent. In order to simplify our analysis, we neglect 
the scale-dependence of cos A; thus, n\ = riz ^ n2 and cos A = sign{B) As/ A. The power 
spectra are written as A|,(A;) = {k^/27r'^){n'^J) = {A^ + A^)^;^"! = A^k^'-'^-^, and A|(A;) = 
{ky2n'^){S^J) = = ^2j2^^n,,„-i Yi&ye defined Uad-l = 2ni and 71^,0-1 = 2^2 to 

coincide with the standard notation for the scalar spectral index. The "isocurvature fraction" 
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defined by fiso = B/A determines the relative amplitude of <S to The cross-correlation 
spectrum is then written as Aj^g{k) — cos A->/A^(A;)A|(A;) = fiso cos A k^'^o,d+niso)/'^-i _ 

The temperature and polarization anisotropics are given by these power spectra: 
oc 

Ci'" oc 

(jcorr ^ 

and the total anisotropy is Cf°* = Cf^ + + 2Cf°''''. Here, gi{k) is the radiation transfer 
function appropriate to adiabatic or isociu^rature perturbations of cither temperature or 
polarization anisotropics. Note that the quantities Uad, niso, and fiso are defined at a specific 
wavenumber ko, which we take to be ko = 0.05 Mpc~^ in the MCMC To translate the 
constraint on fiso to any other wavenumber, one uses 

(I. \ (niso-nad)/'^ 
■ (37) 

We can restrict fiso > without loss of generality. Since we can remove A by normalizing 
to the overall amphtude of fluctuations in the WMAP data, we are left with 4 parameters, 
^ad, niso, fiso, and cos A. We neglect the contribution of tensor modes, as the addition of 
tensors goes in the opposite direction in terms of explaining the low amplitude of the low-/ 
TT power spectrum. We also neglect the scale-dependence of Uad and riiso, as they are not 
well constrained by our data sets. 

We fit to the WMAPext -f-2dFGRS and WMAPext+2dFGRS+Lyman a data sets with 
the 11 parameter model {flth"^, ^mh'^, h, r, Uad, ^iso, fiso, cos A, A, (3, ap). The results of the 
fit for the double infiation model parameters are shown in Table 4. Figure 6 shows the cu- 
mulative distribution of fiso- The best-fit non- primordial cosmological parameter constraints 
are very similar to the single field case. 



/dk j k 
~k\h 



f (f ) [gf{k)gr{k)\ , 



(34) 
(35) 
(36) 
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Table 4. Cosmological Parameters: Adiabatic + Isocurvature Model 



Parameter 


WMAPext+2dFGRS 


WMAPext±2dFGRS±Lyman a 


fi,o{ko = 0.05 Mpc-') 


< 0.32^ 


< 0.33^ 


riad 


0.97 ±0.03 


0.95 ± 0.03 




1 26+°-^^ 

-■-•^"-0.57 


-■-•^^-0.56 


COS A 


'"-0.14 


-n 76+0-18 


A{ko = 0.05 Mpc"^) 


0.82 ±0.10 


0.78 ±0.08 




0.023 ±0.001 


0.023 ±0.001 




0.133 ±0.007 


0.131 ±0.006 


h 


0.072 ± 0.04 


0.072 ± 0.04 


T 


0.16 ±0.06 


0.14 ±0.06 


0-8 


0.84 ±0.06 


0.81 ± 0.04 



^The constraint on the isocurvature fraction, /j^o, is a 95% upper limit. 
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While the fit tries to reduce the large-scale anisotropy with an admixture of correlated 
isocurvature modes as expected (note that cos A < corresponds to IZrad and Srad having 
the same sign, from the definition of initial conditions in the CMBFAST code), this only leads 
to a small reduction in amplitude at the quadrupole. Table 5 compares the goodness-of-fit 
for this model along with the maximum likelihood models for the ACDM and single field 
infiation cases. Because xlff/^^^ ^ot improved by the addition of three new parameters and 
considerable physical complexity, we conclude that the data do not require this model. This 
implies that the initial conditions are consistent with being fully adiabatic. 
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Table 5. Goodness-of-Fit Comparison for Adiabatic/Isocurvature Model 



Model 




ACDM 


1468/1381 


Single field inflation 


1464/1379 


Adiabatic / Isocurvature 


1468/1378 



'^These x^yj values are for the 
WMAPext+2dFGRS data set. Here 
we do not give x^jj for the Lyman 
a data, as the covariance between the 
data points is not known (Verde et al. 
2003). 
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5. 



SMOOTHNESS OF THE INFLATON POTENTIAL 



Spergel et al. (2003) point out that there are several sharp features in the WMAP TT 
angular power spectrum that contribute to the leduced-xlff for the best-fit model being 
~ 1.09. The large xlff may result from neglecting 0.5-1% contributions to the WMAP TT 
power spectrum covariance matrix; for example, gravitational lensing of the CMB, beam 
asymmetry, and non-Gaussianity in noise maps. When included, these effects will likely 
improve the reduced-Xg^j of the best-fit ACDM model. At the moment we cannot attach any 
astrophysical reality to these features. Similar features appear in Monte Carlo simulations. 

While we do not claim these glitches are cosmologically significant, it is intriguing to 
consider what they might imply if they turn out to be significant after further scrutiny. 

In this section we investigate whether the reduced-Xg^j is improved by trying to fit 
one or more of these "ghtches" with a feature in the infiationary potential. Adams et al. 
(1997) show that a class of models derived from supergravity theories naturally gives rise to 
infiaton potentials with a large number of sudden downward steps. Each step corresponds 
to a symmetry-breaking phase transition in a field coupled to the infiaton, since the mass 
changes suddenly when each transition occurs. If infiation occurred in the manner suggested 
by these authors, a spectral feature is expected every 10-15 e-folds. Therefore, one of these 
features may be visible in the CMB or large-scale structure spectra. 

We use the formalism adopted by Adams et al. (2001), and model the step by the 
potential 



where is the infiaton field, and the potential has a step starting at (ps with amplitude and 
gradient determined by c and d respectively. In physically realistic models, the presence of 
the step does not interrupt infiation, but affects density perturbations by introducing scale- 
dependent oscillations. Adams et al. (2001) describe the phenomenology of these models: the 
sharper the step, the larger the amplitude and longevity of the "ringing." For our calculations 
of the power spectrum in these models, we numerically integrate the Klein-Gordon equation 
using the prescription of Adams et al. (2001). 

We also phenomenologically model a dip in the infiaton potential using a toy model of 
a Gaussian dip centered at 0^ with height c and width d: 



We fix the non-primordial cosmological parameters at the maximum likelihood values 
for the ACDM model fitted to the WMAPe^t data, [Qbh^ = 0.022, Q^/i^ ^ o.l3, r = 0.11, 




(38) 




(39) 
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A — 0.74, h — 0.72]. We then run simulated annealing codes for only the three parameters: 
(j)s, c, and d, for each potential, fitting to the WMAP TT and TE data only. For this section, 
since this model predicts sharp features in the angular power spectrum, we had to modify 
the standard CMBFAST splining resolution, splining at A/ = 1 for 2 < / < 50 and A/ = 5 
for / > 50. 

The best-fit parameters found for each potential are given in Table 6, along with the 
xlff for the WMAP TT and TE data. Figure 7 shows these models plotted along with 
the WMAP TT data. The best-fit models predict features in the TE spectrum at specific 
multipoles, which arc well below detection, given the current uncertainties. The step model 
differs from the ACDM model by ^xlff = 10, the dip model by ^xlfj = 6. We are not 
claiming that these are the best possible models in this parameter space, only that these are 
the best-fit models found in 8 simulated anneahng runs. Note that the models with features 
were not allowed the freedom to improve the fit by adjusting the cosmological parameters. 
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Table 6. Best-Fit Models with Potential Features'" 



Model 


0. (Mpi) 


c 


d (Mpi) 


WMAP xlff/iy 


Step 


15.5379 


0.00091 


0.01418 


1422/1339 


Dip 


15.51757 


0.00041 


0.00847 


1426/1339 


ACDM 


N/A 


N/A 


N/A 


1432/1342 



^We give as many significant figures as are needed in order to reproduce our results. 
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A very small fractional change in the inflaton potential amplitude, c ~ 0.1%, is sufficient 
to cause sharp features in the angular power spectrum. Models with much larger c would 
have dramatic effects that are not seen in the WMAP angular power spectrum. 

These models also predict sharp features in the large-scale structure power spectrum. 
Figure 8 shows the matter power spectra for the best- fit step/dip models. Forthcoming large- 
scale structure surveys may look for the presence of such features, and test the viability of 
these models. 



6. CONCLUSIONS 

WMAP has made six key observations that are of importance in constraining infiationary 
models. 

(a) The universe is consistent with being fiat (Spergel et al. 2003). 

(b) The primordial fiuctuations are described by random Gaussian fields (Komatsu et al. 
2003). 

(c) We have shown that the WMAP detection of an anti-correlation between CMB tem- 
perature and polarization fluctuations at ^ > 2° is a distinctive signature of adiabatic 
fiuctuations on superhorizon scales at the epoch of decoupling. This detection agrees 
with a fundamental prediction of the infiationary paradigm. 

(d) In combination with complementary CMB data (the CBI and the ACBAR data), the 
2dFGRS large-scale structure data, and Lyman a forest data, WMAP data constrain 
the primordial scalar and tensor power spectra predicted by single-field inflationary 
models. For the scalar modes, the mean and the 68% error level of the 1-d marginalized 
likelihood for the power spectrum slope and the running of the spectral index are, 
respectively, Usih = 0.002 Mpc~^) = 1.13 ± 0.08 and dus/dlnk = -0.055l^:^i. This 
value is in agreement with dus/dlak — — 0.0311q;qi8 of Spergel et al. (2003), which 
was obtained for a ACDM model with no tensors and a running spectral index. The 
data suggest at the 2-o" level, but do not require that, the scalar spectral index runs 
from Us > 1 on large scales to < 1 on small scales. If true, the third derivative of 
the inflaton potential would be important in describing its dynamics. 

(e) The WMAPext+2dFGRS constraints on Ug, dus/dlnk, and r put limits on single-field 
infiationary models that give rise to a large tensor contribution and a red {ug < 1) tilt. 



-35- 



A minimally-coupled \(f)^ model lies more than 3-a away from the maximum likelihood 
point. The contribution to the Ax^ between the two points from WMAP alone is 14. 

(f) We test two-field inflationary models with an admixture of adiabatic and CDM isocur- 
vature components. The data do not justify adding the additional parameters needed 
for this model, and the initial conditions are consistent with being purely adiabatic. 

WMAP both confirms the basic tenets of the inflationary paradigm and begins to quanti- 
tatively test inflationary models. However, we cannot yet distinguish between broad classes 
of inflationary theories which have different physical motivations. In order to go beyond 
model building and learn something about the physics of the early universe, it is important 
to be able to make such distinctions at high significance. To accomplish this, one require- 
ment is a better measurement of the fiuctuations at high /, and a better measurement of r, 
in order to break the degeneracy between Ug and r. 

We note that an exact scale-invariant spectrum (n^ = 1 and dns/dlnk = 0) is not yet 
excluded at more than 2a level. Excluding this point would have profound implications in 
support of infiation, as physical single field infiationary models predict non-zero deviation 
from exact scale-invariance. 

We conclude by showing the tensor temperature and polarization power spectra for the 
maximum likelihood single-field infiation model for the WMAPext+ 2 dFGRS+ Lyman a data 
set, which has tensor/scalar ratio r — 0.42 (Figure 9). The detection and measurement of 
the gravity-wave power spectrum would provide the next important key test of infiation. 
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at NASA Headquarters and by the hard and capable work of scores of scientists, engineers, 
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and figures. We thank Uros Seljak for his help with modifications to CMBFAST. HVP 
acknowledges the support of a Dodds Fellowship granted by Princeton University. LV is 
supported by NASA through Chandra Fellowship PF2-30022 issued by the Chandra X-ray 
Observatory center, which is operated by the Smithsonian Astrophysical Observatory for 
and on behalf of NASA under contract NAS8-39073. We thank Martin Kunz for providing 
the causal seed simulation results for Figure 1 and Will Kinney for useful discussions about 
Monte Carlo simulations of ffow equations. 
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A. INFLATIONARY FLOW EQUATIONS 

We begin by describing the hierarchy of inflationary flow equations described by the 
generahzed "Hubble Slow Roll" (HSR) parameters. In the Hamilton-, J acobi formulation of 
inflationary dynamics, one expresses the Hubble parameter directly as a function of the field 
(p rather than a function of time, H = H{(f)), under the assumption that is monotonic in 
time. Then the equations of motion for the field and background are given by: 



4> = -2Ml,H\(t>), 
1 



2M2 



V{<t>). 



(Al) 
(A2) 



Here, prime denotes derivatives with respect to 0. Equation (A2), referred to as the 
Hamilton- Jacobi Equation, allows us to consider inflation in terms of rather than 

V{(j)). The former, being a geometric quantity, describes inflation more naturally. Given 
H{(f)), equation (A2) immediately gives V{(f)), and one obtains H(t) by using equation (Al) 
to convert between H' and H. This can then be integrated to give a{t) if desired, since 
H{t) = d/a. Rewriting equation (A2) as 



1 - -e« 
3 



(A3) 



we obtain 



[V{cj>) - 02 



so that the condition for inflation (a/a) > is simply given by e^^ < 1. 

Thus, one can deflne a set of HSR parameters in analogy to the PSR parameters of 
§ 3.2.2, though there is no assumption of slow-roll implicit in this deflnition: 



Vh = 



2Ml 



(2Mpi)' 



EM 

H{ct>) 
H'{4>)H"'{4>) 



(A4) 
(A5) 
(A6) 
(AT) 
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We need one more ingredient; the number of e-folds before the end of inflation, N is defined 



N= Hdt = 

Jt Jd) 



(A8) 



where and 0e are the time and field value at the end of inflation, and N increases the 
earlier one goes back in time (t > =^ dN < 0). The derivative with respect to N is 
therefore, 

W A/f . 

-V~e^. (A9) 



_d_ 
dN 



2 " 'd(f) 

Then, an infinite hierarchy of inflationary "flow" equations can be defined by differentiating 
equations (A4)-(A7) with respect to N: 



den 
dN 

dN 



= [(£ - l)r,, - {'K) + (£ > 0) 



(AlO) 

(All) 



The definition of the scalar and tensor power spectra are: 



- 



H 



2-K J 



k=aH 



8 



(A12) 

(A13) 



Since derivatives with respect to wavenumber k can be expressed with respect to N as: 

d 



d 
dN 



din A;' 



(A14) 



the observables are given in terms of the HSR parameters to second order as (Stewart & 
Lyth 1993; Liddle et al. 1994), 



r = 16e« [1 + 2C(e^ - ?7^)] 



Us 



1 = (2r^, 



4e, 



dus 
dink 



l--(3-5C)e, 

dus 
'dN' 



(A15) 

(3_5C)e^ + i(3-C)e« (A16) 



(A17) 



where C = 4(ln2 + 7) — 5 and 7 ~ 0.577 is Euler's constant. Note that, as pointed out in 
Kinney (2002b), there is a typographical error in defining C in Liddle et al. (1994) that was 
inherited by Kinney (2002a). We have used the correct value from Stewart & Lyth (1993). 
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Finally, the PSR parameters are given in terms of the HSR parameters to first order in 
slow roll as: 

tH = ey (A18) 
tIh = Vv - ev (A19) 
e« - ^y-3eyriy + 3el. (A20) 
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Fig. 1. — Temperature-Polarization angular power spectrum. The large-angle TE power 
spectrum predicted in primordial adiabatic models (solid), primordial isocurvature models 
(dashed), and in causal scahng seed models (dotted). The WMAP TE data (Kogut et al. 
2003) is shown for comparison, in bins of Al — 10. 
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Fig. 2.— This figure shows as a function of k for the WMAP (left), WMAPext+2dFGRS 
(middle) and WMAPext+2dFGRS+Lyman a (right) data sets. The mean (solid line) and 
the 68% (shaded area) and 95% (dashed lines) intervals are shown. The scales probed by 
WMAP, 2dFGRS and Lyman a are indicated on the figure. 




Fig. 3. — This set of figures shows part of the parameter space spanned by viable slow roll 
inflation models, with the WMAP 68% confidence region shown in dark blue and the 95% 
confidence region shown in light blue. 
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Fig. 4. — This set of figures compares the fits from the WMAP (top row), 
WMAPext+2dFGRS (middle row) and WMAPext+2dFGRS+Lya data (bottom row) to 
the predictions of specific classes of physically motivated inflation models. The color coding 
shows model classes referred to in the text: (A) red, (B) green, (C) magenta, (D) black. 
The dark and light blue regions are the joint 1-cr and 2-0" regions for the specified data sets 
(contrast this with the 1-d marginalized 1-cr errors given in Table 1). We show only Monte 
Carlo models that are consistent with all three 2-0" regions in each data set. This figure does 
not imply that the models not plotted are ruled out. 



- 46 - 




TO M 1.4- ij! 1.0 J.J ij i.t- 030 o;a Oj« o.e u iji 

Iktll.R-ajK£M(!t'^ <kdfc4'i>m'Uii«'' InWtli'H^ldbr- 




!l I:^ 1.4 14 O-S 1.(1 I. Id 1.4 1.E DjQ tl,;t .11.4. Cut OA 1J): Jj$ 

■! kii[UXt£M]ic' i),d1i!30.<]<J'i! h^' linaurriic^iitiiDr 

Fig. 5. — This set of figures compares the fits from the WMAPext + 2 dFGRS+Ly a data to 
the predictions of all four classes of inflation models. The top row is Class A [red dots] . The 
second row is Class B [green dots] . The third row is Class C [magenta dots] . The bottom row 
is Class D [black dots]. The dark and light blue regions are the joint l-a and 2-0" regions for 
the WMAPext + 2 dFGRS+Ly a data. We show only Monte Carlo models that are consistent 
with 2-0" regions in all panels. This figure does not imply that the models not plotted are 
ruled out. 
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Fig. 6. — The cumulative distribution of the isocurvature fraction, fiso, for the 
WMAPext+2dFGRS+Lyman a data set. 
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Fig. 7. — Best-fit models (solid) with a step (left) and a dip (right) in the infiaton potential, 
with the WMAP TT data. The best-fit ACDM model to WMAPext data is shown (dotted) 
for comparison. 




Fig. 8. — The large-scale structure power spectra for the best-fit potential step (left) and 
dip (right) models. 
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Fig. 9. — The tensor power spectrum for the maximum hkehhood model from a fit to 
WMAPext+2dFGRS data sets. The plot shows the TT (solid), EE (dots), BB (short dashes) 
and the absolute value of TE negative (dots and dashes) and positive (long dashes) tensor 
spectra. 



